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CN . Abstract. After a brief discussion of the meaning of the potential in quantum mechanics, 

we shall examine the results of the Yukawa model (scalar meson exchange) for the nucleon- 
nucleon interaction in three different dynamical frameworks: the non-relativistic dynamics 
of the Schrodinger equation, the relativistic quantum mechanics of the Bethe-Salpeter and 
Light-Front equations and the lattice solution of the Quantum Field Theory, obtained in the 
quenched approximation. 



" 1 Introduction 

Since Newton's time, an interaction between two particles is understood as something that pre- 
1 vents their relative motion to be rectilinear and uniform. In non-relativistic Quantum Mechanics 
this is realized by any operator V, traditionally called potential, that disturbs a plane wave |A;), 
i.e. there is no any constant A such that 



(H + V)\k) = \\ k) 



while 



E being the energy of the state \k). 

From where does the potential V come from? Sometimes it is taken from classical mechanics, 
like in the Coulomb case, but it can be picked out of a hat as well, provided one respects some 
space-time (translation, rotation, P, T) or internal (C, isospin) symmetries. 

This way to built an interaction is perfectly legal and can be even extended in order to satisfy 
the requirements of relativistic invariance. The proper method of constructing the Poincare 
algebra with a given V was formulated in the series of works by Bakamjiam and Thomas PQ, 
extended later by Keister and Polyzou [2]. 

From this point of view, particles interact because they go into an inhomogeneous region 
of the space: if V was a constant, nothing would happen. This "ex nihilo" approach, though 
leading to some remarkable success, has not been very fertile when describing the physical 
phenomena, specially in the subatomic world. Here is the crucible where deep changes in the 
states of matter occur, one of the deepest being the non-conservation of the particle number, 
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like in the annihilation of matter into light, e + + e~ —7-7 + 7, or > conversely, the transform of 
"speed" into matter-antimatter pairs like in p + p -^-p + p + p + p, both processes which are 
nowadays customary in the laboratories. The potential approach of the non-relativistic quantum 
mechanics with a fixed number of particles is here of a little help. Even worst, it does not provide 
any "way of thinking" about the natural processes. 

A more interesting approach is provided by the Quantum Field Theory (QFT), now custom- 
ary in the high energy physics world. In QFT, the interaction is a consequence of the exchange 
of a bosonic mediator field. This approach has been successfully applied since the development 
of quantum electrodynamics to every piece of the Standard Model of particle physics. In the 
Lagrangian formulation, the simplest case is provided by the Yukawa model j3]- The interaction 
between fermionic matter fields \P is mediated by a scalar field <f> and written in terms of the 
Lagrangian as: 

C(x) = C D [W] + C KG [$\ + dntW, $] 

where: 

Cint(x) = gV(x) <P(x) 9{x) 
and Cd, £kg are respectivley the Dirac and Klein-Gordon free Lagrangians. 

Within this framework, a particle on a free state \k\) emits a quanta \q) which is absorbed 
by a particle on a state \k-zj. Their initial states are modified in the process and results into new 
ones \k'i) and \k' 2 ) : they have "interacted". 




Figure 1. Lowest order exchange graph. 



The basic element of this "lego" is the annihilation-creation-creation in the interaction vertex, 
driven by a strength constant go : 

9 a{, b\ a K -> g a^e"^* 1 b\e^ a kl e ik ^ -> g V(xi)$(xi)V( Xl ) 

and Quantum Field Theory tells us how to associate to the process displayed in Fig. [T]a proba- 
bility amplitude A(ki, &2 — > k[, k' 2 ). 

The relation between the QFT and the potential approach is made by identifying V to the 
amplitude of the lowest order "exchange" graph, which according to Feynman rules reads 

V ^ ^ A{klM k '^ k>2) = ^faWJ u{k 2 )u{k' 2 ) . (2) 

The potential V(r) used in non-relativistic quantum mechanics is obtained by a Fourier 
transform of the amplitude A after some - quite crude - simplifications: 



u(ki)u(k'i) = u(k2)u(k 2 ) = 2m 

q = (0,g) = (0,fci-fc 2 ) 
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One thus obtains 



V(k 1 ,k 2 ,k[,k 2 ) = V(q) 



9 



Li 2 +q 2 



and, by a three-dimensional Fourier transform, the usual Yukawa potential 

g 2 e~^ r 
47T r 



V{r) 



(3) 



is obtained. This is the procedure leading to the usual 1/r Coulomb potential starting from 
QED. The same that brought Yukawa to formulate the first theory of strong interactions that 
deserved him a Noble prize. 

Some remarks are in order: 

(i) If the interaction was limited to only one exchange of Fig. [TJ there would never exist a 
bound state. These states are indeed associated to poles in the scattering amplitude, and the 
Born term (one exchange term) has no one. An infinite sum of exchanges is needed to generate 
these singularities and this task is ensured by the dynamical equations 



+ 



+ 



+ 



+ 





(3a) (3b) (3c) 

Figure 2. Perturbative expansion of the scattering amplitude. 



(ii) The ladder sum, used in all nuclear models, accounts only for a small, though infinite, 
part of the full interaction given by £jnt- At any order in the coupling constant, there are several 
contributions which are ignored. Figure [2] illustrates the perturbative series of the scattering 
amplitude: diagrams (2a), (2c) and (2d) are of the same order than (2b) but are not included in 
Eq. (j4]). It is even worse at higher orders, since the fraction of dropped diagrams increases very 
fast. 

Even assuming that the effect of some of the neglected diagrams is incorporated in the 
renomalized quantities - e.g. (2a) in the mass and (2c) in the coupling constant - there remains 
an infinity of them. If g is large, what one neglects is not negligible and one can expect that the 
physics with V could seriously depart from the ones contained in the underlying £; nt . 
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We present in what follows the results obtained when the very same interaction - the Yukawa 
model - is considered in three different dynamical approaches. They will be limited to the 
J n = + bound states and corresponding low energy scattering parameters. 

Apart from being at the origin of the theory of nuclear forces, this model, the simplest 
meson- fermion interaction Lagrangian, has several advantages: (i) in the non-relativistic limit it 
gives the same result for the two- fermion and the two-boson system, (ii) when inserted in the 
relativistic equations - at least the ones considered here - it does not require any regularization 
procedure to be integrated and (Hi) it is a renormalizable quantum field theory. 

Section 2 is devoted to the non-relativistic results. They are widely used and constitute the 
reference ground of most of the nuclear and atomic physics calculations. In section 3 we will 
consider the results of two relevant, relativistic equations: the Bethe-Salpeter [H [5] and Light- 
Front Dynamics [6]. Section 4 will contain the quantum field results obtained using the lattice 
techniques in the quenched approximation. 



2 Non-relativistic results 



We consider in this section the non-relativistic system of two particles with equal mass m, 
interacting by a Yukawa potential (|3|) of strength g and range parameter ji. 

Although the problem depends on three parameters (m, g, (i) it can be shown that the binding 
energy (B) and the scattering length (ao) are given by 

B = m (^) 2 e(G) (5) 

a = -A(G) (6) 
A» 

where e(G) and A(G) are respectively the binding energy and scattering length of the dimen- 
sionless S-wave Schrodinger equation: 



u"{x) + 



-e + G 



u{x) = (7) 



with a coupling constant G, related to the original parameters (m, g, /i) by 

47T fj, 

The functions e(G) and A(G) are displayed in Figs. [3] and H] respectively. The convention 
used for the scattering length corresponds to 5(k) = —a^k + o(k 2 ). The critical value for the 
appearance of the ground state is Go ~ 1.680. At this value A(G) has a pole and it can be shown 
that for small values of G one has 

A(G) = —G + o(G 2 ) (8) 

which corresponds to the Born approximation. All the physics of the non-relativistic problem 
(S-wave) is contained in these two figures, with the understanding that there exists an infinity 
of similar branches in B{G) - and corresponding poles in A(G) - at increasing values of G 
(Go = 1.680, Gi = 6.445, G2 = 14.34,...) indicating the appearance of an infinite number of 
excited states. 

In summary, for a massive exchange /x, there is a non-zero minimal value of the coupling 
constant go required to have the first bound state. It is given by go 2 (m,(i) = 47rGo^ with 
Go ~ 1.680. This go value decreases linearly with ji and vanishes in the Coulomb limit — > 0) 
but in the nuclear case (/x/m ~ 0.5) is rather large go ~ 3. Once the bound state appears, the 
solutions of Eq. (J7|) exist for any value of the parameter G and one can obtain any value for the 
binding energy B, even exceeding 2m, 
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Figure 3. Ground state binding energy of the di- 
mensionless non-relativistic Yukawa model © as a 
function of the coupling constant G. The appearance 
of the first bound state corresponds to Go = 1.680. 




Figure 4. Scattering length in the same model than 
Fig. [3] The Born approximation is indicated by the 
dashed line. The singularity (vertical line) corre- 
sponds to the appearance of the first bound state. 



3 The Yukawa model in Relativistic Dynamics 

Things become less simple when considering the same model in a relativistic framework. This 
covers a very wide and -to some extent- not well defined domain of theoretical physics aiming to 
incorporate all or part of the relativistic invariance in the dynamical equations. It goes from the 
simple implementation of relativistic kinematics to the full Quantum Field Theory treatment, 
which will be the proper way to incorporate relativity to the quantum world but whose solutions 
are very difficult to obtain beyond the perturbative domain. 

We will consider here two of the many relativistic approaches: the Light Front Dynamics 
(LFD) and the Bethe-Salpeter (BS) equation. Although far from being representative of this 
vast domain they illustrate well the kind of qualitative agreements and quantitative differences 
in the predictions they give. Both are rooted in the Quantum Field Theory but present important 
differences in the way they are constructed as well as in the formal objects they deal with. 

The BS equation deals with a field theoretical object corresponding to the following ampli- 
tude [5]: 

<P( Xl ,x 2 , P) = (0\T {W( Xl )§(x 2 )} \P) (9) 

where P 2 = M 2 is the total mass of the two-body system. 

It is usually written in momentum space <P(k,P), obtained after taking into account trans- 
lational invariance and performing a Fourier transform with respect to the relative coordinate 

X = X\ — X2- 

<P(x 1 ,X2,P) = *(x,P) e- iP ^ + ^' 2 

#(x,P) = J J0p*(k,P) e~^. 

The BS amplitude for a two-fermions system is a 4 x 4 matrix in spinor space which can be 
expanded in a basis of independent Dirac structures S( c > 

n c 

<2> = X> S (c) S® € {1,7^75,7^75,^} (10) 
c=l 
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The number of scalar components n c depends on the quantum number of the state. For J = + 
there are four of them which can be chosen as: 

P) = Mk, PH + Mk, P)^JB + Mk, P) 

where k = 7^ and = §[7^,7^]. 

The dynamical equation can be represented in the following graphical form: 



(k-P)P 



k_ 

M 



-/ 5 + Mk,P) 



^fj,uPfj,k v 



M 2 



(11) 




For the case of two equal mass fermions it reads 

S^(h)<P(k,P)S^(k 2 ) = J iK(k,k') A $(k! , P) r 2 (12) 

where 

i P 

Si(h) = j , k 1 = — +k 

ki — m + ie * 

S2{k 2 ) = —~ , k 2 = — -k 

—k 2 — m + ie * 

are the fermion propagators, and K are respectively the vertex functions and the meson 
propagator. They depend on the particular type of meson-fermion coupling. In the case of the 
Yukawa model (scalar meson exchange) they are 

"^■"'- FiFT^ - (13) 

and 

A = r 2 = r 2 = -ig 



The BS equation (|12[) in momentum space is four dimensional. After a partial wave expansion 
it reduces to a set of coupled two-dimensional equations among the different components 4> c of 
the state (|10p. This equation has been solved by several authors for the bound state problem 
both in the Euclidean [8] and Minkoswki metric [£l [10] and for the on-mass shell scattering 
amplitudes [11 ( 112 1 IT5]. For the scattering states, the full (off-shell) solution in Minkowski space 
has been obtained only very recently |14[ [T5| ITS] . 

LFD can be understood as an Hamiltonian formulation of the QFT defined on an space-time 
surface of equation ui-x = a, where w is a light cone vector uJ 2 = [6j. 

The state vector {^(a)) is defined on this plane and the Poincare algebra generators are 
obtained by integrating through this surface the flux of the conserved Noether currents associated 
to a given Lagrangian C. In case of translations, for instance, they are given by 

P»(a) = [ T^(x)5{u-x - a) uj v d 4 x (14) 
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dC 



d v <j) - g^C 



and fulfill 



d a P»{a) = 



Once the generators are obtained, the dynamical equation determining the mass of the system 
M 2 is given by 



P 2 | <F) = M 2 | <P) 



After some algebra, one is led to 



(M 2 - P 2 ) | V) = 2P -lo J Hint^r) exp(-mr) dr \ V) 



(15) 



(16) 



where 7-Li n t{k) denotes the Fourier transform of the hamiltonian density 

««(*) = /««(x)e X p(«.,)A 

The state vector is decomposed into its Fock components with an increasing number of 
particles, which can be schematically written as: 



9) 



[ d A k x . . . d 4 k a <P qi ...d A 

a/3 J 



9a/}(ki, ...,k a ,qi...q p )al l ... a\ a b\ ± . . . b\ | 0) 



(17) 



The components of this expansion W a a are the relativistic counterparts of the usual non- 
relativistic wave functions. They have also a probability interpretation and are smooth func- 
tions of the arguments. We can consider the set ^n al3 = {^afi} as the components of an infinite 
dimensional vector & = (l^q, <p2) ^3> • • •), coupled to each other via the interaction operator Hint- 
Equation (fTBj) is thus an infinite system of coupled channels. 



(M 2 - Pi 



f...\ 

\...J 



2P -lj J Hint(ur)dT 



f...\ 

#2 



If we restrict ourselves to the two- = {^20}) and three-body (1^3 = {$21}) wave functions, 
we obtain a system of two coupled equations for W 2 and ^3 which constitutes the ladder ap- 
proximation. By expressing 1F3 in terms of one gets an integral equation for ^ with an 
energy-dependent kernel. 

The LF equation for the two-fermion system is three-dimensional and takes the form 



[M 2 — 4(fc 2 + m 2 )] <P(k,n) 



^ jK(k,k',n,M 2 )<P(k',n)^- 



(18) 



where K(k, k ', n, M 2 ) is the interaction kernel and n is unit vector, the spacial part of the light 
cone one ui = (u>o,u}on). For the Yukawa model it reads 



K(k,k',n,M 2 



Am 2 (Q 2 + /i 2 ) 



u{k 2 )u a > 2 (k' 2 ) u(/ci)u CT / (k[ 



(19) 
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with 

Q 2 = (k - k'f - (n-fe)(n-fe') {£k ' ~ £k? + U, +e\- hA 

£k'£k V 1 J 

The J = + wave function has the form 

<P(k, n) = u(k 2 )(/)Ujl(k 1 ), (21) 
where eft is expanded in terms of spin structures Si 

<t> = fiSi+f 2 S 2 , (22) 

61 = o /n 75, 5 2 = —7= — — 2" 75 23 

2V2efc 2V2m | k x n | \ oj-p e% J 

and /j are scalar components depending on (k, k n) . 

By inserting the expansion (|21[) in (|18p it results - like for the BS case - in a system of 
two-dimensional integral equations coupling the different components of the wave function fa 
despite the fact that the LF equations is only three-dimensional. This is due to the existence of 
an additional vector n in the theory. Notice however that the number of components in LFD 
is 2, half the number in the BS case. The explicit equations and the numerical solutions of the 
LFD equation for the Yukawa model have been presented in[I7lHglHgil5Dll2Ill221l251l23|. 



nk nk 



(20) 




Figure 5. Binding energy of the J — + state in Yukawa model as a function of the coupling constant g 2 given by 
the Light-Front (LF), Bethe-Salpeter (BS) and non-relativistic Schrodinger equation. Vertical dotted line denotes 
the critical coupling constant g c = 6.84 for LF and g c = 2tt for BS. 



The BS and LFD binding energy of the J = + state as a function of the coupling constant is 
displayed in Fig.[5]for the parameters m = 1 and fi = 0.5 and is compared to the non-relativistic 
results given by ©. 

It is worth noticing the existence, for both relativistic equations, of a critical coupling con- 
stant g c . Above that value, g > g c , the systems "collapses", i.e. its spectrum is unbounded, and 
vertex form factors are required to solve the corresponding dynamical equations. Indeed, when 
g — > g c from below, the value of M 2 vanishes, becomes negative and tends smoothly to —00. 
This happens in a rather narrow domain where g is very close to g c and which is not very dis- 
tinguished in Fig. [5l The physical meaning of BS solution is lost already at g < g c though very 
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close to, when M 2 becomes negative. This result was first found in [191 [20"1 EI] in the framework 
of LFD and also obtained for the BS equations [10] using the methods developed in [20 . 

If the very existence of this critical coupling constant is common to both relativistic ap- 
proaches, their precise numerical value, independent of m and //, however depends on the par- 
ticular dynamics: one has g c = 2ir for BS and g c = 6.84 for LFD. This difference is due to the 
different treatment of the intermediate states in the ladder kernel in these two approaches: while 
the ladder BS equation incorporate effectively the so-called stretch-box diagrams [25] , they are 
absent in the ladder LFD results. 

Contrarily to the non-relativistic case, the range of the strength parameter g in these rela- 
tivistic equations is limited. These limits are indicated by vertical dashed lines in Fig. [5l As one 
can see the accessible binding energies are the same in both equaztions B £ [0, M] due to the 
vanishing of M 2 near g c . 

The results of both relativistic equations are quite close to each other, specially at moderate 
values of B (see right panel of Fig. [5]) but depart from the non-relativistic ones which, for a 
given value of g, generate always much more attraction. This differences are not of kinematical 
origin, since they exist even in the limit of zero binding energies and increase with the value of 
\x. Strictly speaking the results would coincide only in the limits (i — > (Columb problem) and 

To our knowledge there are no published results for the scattering observables with the 
Yukawa model and the considered equations. They have been however computed for the scalar 
model ((f) 2 x theory) both in LF [261 [27] and in the BS one [281 HH El CGS UM- 

The results described above illustrate well the kind of dispersion one can find when moving 
from a non-relativistic to a relativistic description of the same system, would be the simplest one 
and submitted to the simplest interaction. While the result of the non-relativistic Schrodinger 
equation with a given potential is unique, it is not the case in the relativistic world. The imple- 
mentation of relativity can be done following different approaches, but they give rise to qualita- 
tive results which are common to some of them. For instance: the strong repulsive effects, the 
existence of critical coupling constants or, when solving the three body system, the automatic 
generation of three-body forces |30[ 131] , The reason for such differences is not in kinematics. 
Notice that, as a consequence of the inequality 



\lfp- + m 2 — m < 0, 

VF 2m 

the relativistic kinematical corrections are always attractive, while the results of figO and similar 
one for the scalar theories, shows rather a strong repulsion. The origin of this new behavior is 
thus dynamical and lies in the interaction kernel as it can be seen by computing the zero energy 
cross sections (see Fig. 1 from [15]) . 

One of the more consistent approaches to relativistic ab initio nuclear physics is the one 
developped by Gross and Stadler using the spectator equation [36J. The philosophy is quite 
close to the BS and LF equation: using this relativistic equation and OBE kernels, these authors 
obtain a very good fit to pn data with a relatively small number of parameters and reproduce the 
experimental triton binding energy without explicitly adding 3-body forces. We would however 
remark, that there exists other relativistic approaches which substantially differ from the ones 
described above. Of particular interest is the approach developped by H. Kamada, W. Gloeckle, 
H. Witala, J. Golak, Ch. Elster, W. Polyzou and coworkers. The starting point is a potential 
which in the non relativistic dynamics provides a satisfactory description of NN data. Using 
the Bakamjiam-Thomas construction [H [32], a new relativistic potential is obtained in such a 
way that once inserted in a relativistic Lipmann-Schwinger equation it produces the same phase 
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shifts than the non relativistic ones [33]. The parameters of this new potential are not readjsuted: 
it is an implicit function of the preceding ones and contains no new parameters. In this scheme 
there are, by construction, no any two-body relativistic effects. The real difference between 
relativistic and non relativistic dynamics appears only when going to the three-body problems. 
This approach has been succesfully applied to the few-nucleon problem |34[ [35] [371 138| [39] . 

4 The Yukawa model in the Lattice 

The very large effects we found when including the cross ladder kernel in the BS and LF equations 
|40j , as well as the pionneer work of [41] computing the full cross ladder sum in the scalar theories 
motivated a work to evaluate the full QFT content of the Yukawa model. 

A series of papers [42 [ |43| [44] has been devoted to this project with the aim of obtaining 
the B(g 2 ) dependence of the Yukawa model as well as some low energy parameters. They are 
summarized in [45] 

To this aim we have used the standard lattice techniques, developed in the framework of 
QCD, and the Lagrangian density: 



C = C D {* 1 *) + C KG {$) + C I {y,y,$) , (24) 



with, in Euclidean space, 



C D {W,W) = Wr(d^ + m )Wr , (25) 

£ KG ($) = i(d M <2W + ^ 2 ) , (26) 

= go&m (27) 

where & denotes respectively the fermion field and <fr the exchanged meson field responsible for 
the interaction. 

The fermion field is supposed to describe a nucleon (N) and the meson field a - more or 
less fictitious - scalar particle (a) responsible for the attractive part of the NN potentials. 
The Lagrangian depends on three parameters: the fermion tuq and meson /io masses and a 
dimensionless coupling constants go. 

The theory is solved in a discretized space-time Euclidean lattice of volume V = L 3 x T 
and lattice spacing a using the Feynman path integral formalism. In this approach, the vacuum 
expectation values of an arbitrary operator O is given by the integral: 

(0(Wr,¥,$)) = | J[d&][d&][d$]0(?,¥,$) e - s *[*.*.*] , (28) 

where, according to (|24p . the discretized Euclidean action Se can be written in the form 

Se = « 4 C = Sd + Skg + Si 

X 

and plays the role of a probability distribution in a Monte Carlo simulation. 

The fermionic part (Sd+Si) is written as a bilinear form in the dimensionless fermion fields 

Sd + S T = J2i> x D xy % (29) 



where 



xy 



D X y = 5 X< y - ft^J [(1 - Jjj) 5 X) y-n + (l + 7^ ) 5 Xi y +(i ] + § L 4> (30) 
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is the Dirac- Yukawa operator, 

1 



+ 2amQ 



(31) 



the hopping parameter and gi = 2k$o the lattice coupling constant. In terms of the dimensionless 
meson field <b = a<l>, the discrete Klein-Gordon action reads: 



(8 + aVo) 0x-2j^' 



(32) 



The integral over the fermion fields in (|28p is performed by algebraic methods. The keystone 
in a lattice simulation is the fermion propagator S(x,y), corresponding to 0(^,ip,(f)) = i/j x ifty 
After performing the fermionic integration, one is left with 

S(x,y) = (^ y ) = iy [#] D~ y det[D(<p)] e~ s ^ , (33) 

This implies the evaluation of a determinant and inverse the Dirac operator that, even for 
moderate lattices V ~ 24 , has a dimension of ~ 10 6 . Moreover, if a Monte Carlo simulation 
is to be done using Eq. (|33p . the probability distribution for meson configurations is given by 
e -S , M(</')-iog(dct(D))^ w j ia t means evaluating a large determinant in every Monte Carlo step . This 
can be avoided by the use of Hybrid Monte Carlo techniques that nevertheless are the main 
source of time spent in the simulation. This task is considerably simplified in the "quenched" 
approximation that, from the computational point of view consists in setting det(D) independent 
of the meson field in the fermionic integral. 

From a physical point of view, the quenched approximation avoids the possibility for a meson 
to create a virtual nucleon-antinucleon pair <P — > (see Fig. [6]). Due to the heaviness of the 
nucleon with respect to the exchanged meson this approximation is fully justified in low energy 
nuclear physics and implicitly assumed in all the potential models. Under this hypothesis the 



Figure 6. The quenched approximation neglects the possibility for a meson $ to create a virtual fermion- 
antifermion pair 

generation of meson- field configurations according to the probability distribution e~ Su ^ is also 
greatly simplified and the path-integral sum over the mesonic fields in [33] can be accurately 
computed, to compute D~ y [<j>] for an statistical ensemble of meson field configurations. 

1 N 

i=l 

Due to translational invariance one is left in practice to compute S(x, 0) = D~ [</)], that is to 
solve the linear system: 



D zx ((p)S x = S z0 



(35) 
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One can obtain in this way the renomalized fermion mass am as well as the mass of the 
two-body interacting particles C1M2, both in lattice units. The binding energy - in constituent 
mass units - is then given by B/m = ((1M2 — 2am) / am. 

In Fig. [8]we show this binding energy as a function of the lattice size Lafi for a given set of 
parameters. The dotted line is a fit obtained with a 1/L 3 dependence. As it can be seen in this 
figure, the binding tends to zero in the infinite volume limit. This indicates that this two-fermion 
system has no bound state for this particular set of parameters. The binding energy results from 
setting two interacting particle in a box with periodic boundary conditions but contrary to a 
real bound state, this one disappears in the limit L — > 00. It turns out that the situation is 
however the same for the whole range of parameters accessible in the numerical simulations. 

As in the relativistic dynamics, though for a completely different reason, there exist a max- 
imum value of the coupling constant that can be attained within this framework. The reason 
is the existence of zero modes in the Dirac- Yukawa operator, i.e. the appearance of meson field 
configurations such that det[-D(^)] ~ thus leading to an ill-conditioned linear system (|35p 
and the impossibility to compute the fermion propagator. As a practical measure of the "ill— 
conditioness" of D we have considered its "condition number" defined as the ratio between the 
largest to the lowest eigenvalue modulus. The largest is this number the more difficult is to solve 
the linear system. Depending on the method used for that purpose, either the algorithm cannot 
find the solution, or the round-off errors make the solution wrong. 

It was found that such "ill-conditioned configurations" appear in the Yukawa model for 
almost any k when gi > 0.6. In this case the inversion of the Dirac operator becomes in practice 
impossible. For illustrative purposes, we have plotted in Fig. [7] the condition number of D as a 
function of the lattice coupling constant gi for an ensemble of L = 8 configurations at fixed value 
of k. As one can see, the condition number of a given configuration diverges on a discrete set of 
gi values for g^ > 0.6 indicating the practical impossibility to compute the nucleon propagator. 
The precise gi values where this divergence occurs depend on the particular configuration, on 
the values of k and afi and on the lattice size. It turns out however that the situation described 
in Fig. [7] is generic for the quenched Yukawa model. 

1e+12 | , , , , 1 



1e+10 - 




0.2 0.4 0.6 0.8 1 

9L 



Figure 7. Conditioning number as a function of gL for a fixed value of k = 0.11 and V = 8 4 and for 9 different 
meson field configurations. 

It is worth noticing that in the full QFT formulation every configuration is weighted by the 
determinant of the Dirac operator D and therefore the configurations yielding an ill-conditioned 
linear system (|35p . i.e with det(D) ~ 0, do not contribute to the functional integral. In the 
quenched approximation, however, this is no longer true and "ill-conditioned configurations" 
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can be sampled. 

Therefore, the numerical simulations in the quenched Yukawa model are limited to values 
of the lattice coupling constant g^ < 0.6. Using a typical value of k = 0.1, this corresponds 

2 

to g = %h < 3, that is a = f- < 0.7 which is of the same order than the «qcd m the non- 
perturbative region. 

Although assuming that this problem could be associated to the quenched approximation it 
is physically surprising that no any NN bound state could be generated if the NN pair creation 
is not taken into account. 
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Figure 8. Binding energy versus lattice size L, for 
ql = 0.3, K — 0.118, and a/j, = 0.1 averaged over 
4000 samples for L = 16, • • • , 32, 2000 for L = 36, 
and 800 for L = 48. Dotted line corresponds to a 
1/L 3 fit. 




G 

Figure 9. Scattering length vs G for a lattice vol- 
ume La[i = 2.4. Solid line indicates the continuum 
non-relativistic result and the dotted one the Born 
approximation |(HJ. Blue circles represent the NR 
results. 



In absence of bound states in the quenched approximation, we still can access to the NN low 
energy scattering parameters. The scattering observables cannot be obtained in Euclidean time 
in the infinite volume limit [16] but can be extracted from the volume dependent binding energy 
measured on finite lattices, like for instance the one plotted in Fig. [8l The underlying formalism 
was developed by Luscher in [37] who gave a 1/L expansion of the binding energy. In its leading 
order it reads: 

m = ~T~\* (36) 

Using the binding energy values of Fig. [8] and equation (|36p . the NN scattering lengths ao 
can be computed. The result - corresponding to gi = 0.3, k = 0.118, and a/j, = 0.1 - is 
ao/x ~ —0.13 and the dimensionless coupling constant of the non-relativistic model is G = 0.193. 
The corresponding non-relativistic scattering length value, given by Fig.U is Aq = a^fi = —0.214. 

This study has been performed for several values of g^. The dependence of ao on the coupling 
constant G is plotted in Fig. O for a lattice size of Lafi = 2.4 (L = 24, ajjL = 0.1). One can see 
that the lattice results notably departs from the non-relativistic ones (solid line) and are above 
the Born approximation (dashed line). 
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The values of the accessible coupling constants extend beyond the Born regime but are still 
far from the pole behavior corresponding to the appearance of the first bound state displayed 
in Fig. SJ The difference between the lattice and NR results may indicate strong repulsive 
corrections. These kind of corrections were already manifested in the bound state problem when 
solving the same Yukawa model both in Light Front and BS ladder equations. 

5 Conclusion 

We have presented the results of the Yukawa model for the NN system using three different 
dynamical frameworks. In this model - which constitutes the simplest renormalizable quantum 
field theory of a fermion-meson interaction - two fermions of identical mass m interact by 
exchanging a massive scalar particle of mass [i. Results are restricted to the J n = + state. 

In non-relativistic dynamics the model consists in solving the Schrodinger equation with the 
so called - static and local - Yukawa potential ([3]). The results depend on three parameters 
(m,(x,g) but some simple scaling properties make it dependent on a single dimensionless param- 
eter, the coupling constant G = (g 2 /47r)(m//iz). The existence of a first bound state requires the 
coupling constant to be greater than a critical value Go = 1.680. For G > Go, its binding energy 
increases monotonously and can reach any arbitrary value, even greater than 2m. 

The way for implementing relativistic dynamics in the description of the same system is 
not unique. We have considered two relevant relativistic equations: Light-Front and Bethe- 
Salpeter. In both cases, for this particular coupling and state, the corresponding equations 
can be integrated without introducing any regularization form factor. The scaling properties 
of the non-relativistic model model are now lost and both equations exhibit a critical coupling 
constant g c , for which M 2 tends to — oo. For slightly smaller value of g the mass M 2 crosses 
zero: the system "collapses". The precise numerical value of g c depends however on the particular 
dynamical framework: g c = 6.84 for Light-Front equation and g c = 2tt for BS one. This value is 
however large enough to generate several bound states. The corresponding relativistic binding 
energies B(g 2 ) are close to each other but they depart sizeably from the non-relativistic ones 
even for very loosely bound states and are strongly repulsive. 

Relativistic equations are the first step towards a full Quantum Field Theory solution. They 
suffer from two main drawbacks. On one hand, most of the one-boson exchange kernels require 
to be regularized in order to obtain an integrable equation. This is usually done by introducing 
a vertex form factor cutting the high momentum components above some arbitrary value A, 
but thus diluting all the benefit of an approach starting from the first principles, like underlying 
Lagrangian. The Yukawa model (in J = + state) is rather an exception than the generic case 
of one-boson exchange models. On the other hand, the ladder kernel accounts only for a small 
part of the interaction, specially when large values of the coupling constant are involved. We 
have presented the first attempt to incorporate the full dynamical content of the Yukawa model 
by using standard lattice techniques developed in the context of QCD. 

The only approximation in the calculations was to neglect the NN loops (quenched approx- 
imation) as it is the case in all the nuclear models. The lattice results indicate the existence 
of a maximal coupling constant gi which is well below the required value for the existence of 
the first bound state in the heory. For smaller coupling constant the scattering length has been 
computed using Luscher's prescription. Results are in agreement with the repulsive effect found 
by solving the relativistic equations. 

We have shown how the same interaction model between fermions can give rise to very differ- 
ent physical pictures depending on the particular dynamical framework in which it is considered. 
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The differences between a relativistic and a non-relativistic approach are not only quantitative 
but lead to new qualitative behavior. In the particular case of the Yukawa model, which was the 
first of all the nucleon-nucleon models, the full quantum field solution remains still unknown. 
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